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Flat band analogues and flux driven extended electronic states in a class of 
geometrically frustrated fractal networks 
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We demonstrate, by explicit construction, that a single band tight binding Hamiltonian defined 
on a class of deterministic fractals of the b — 3N Sierpinski type can give rise to an infinity of 
dispersionless, fiat-band like states which can be worked out analytically using the scale invariance 
of the underlying lattice. The states are localized over clusters of increasing sizes, displaying the 
existence of a multitude of localization areas. The onset of localization can, in principle, be ‘delayed’ 
in space by an appropriate choice of the energy of the electron. A uniform magnetic field threading 
the elementary plaquettes of the network is shown to destroy this staggered localization and gener¬ 
ate absolutely continuous sub-bands in the energy spectrum of these non-translationally invariant 
networks. 
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I. INTRODUCTION 


The interplay of geometry and quantum interference 
leading to exotic band structure in a class of quasi 
one-dimensional geometrically frustrated lattices (GFL) 
has drawn considerable interest in recent times. The 
GFL have been quite exhaustively studied over the years 
mainly in the context of localized spin systems [l| , using 
the Heisenbergmodel for studying frustrated antiferro¬ 
magnetism 0 JlL or itinerant ferromagnetism within the 
Hubbard model 1,0. The degeneracy induced by the ge¬ 
ometric frustration and non-trivial ground states such as 
the spin liquid states as its consequence [l| have unveiled 
rich physics of the spin systems. 

Parallely, electronic spectra of several simple lattices, 
studied mainly within a single band tight binding for¬ 
malism and only with nearest neighbor hopping, have 
displayed exotic features such as the appearance of dis¬ 
persionless, flat bands (FB) [I3-0- Such frustrated 
hopping models are important as they offer prospects of 
strong interaction physics such as the fractional quantum 
Hall effect. Graphene anti-dot lattices [ll| and post¬ 
graphene materials such as the phosphorene [ll| have 
been proposed as systems exhibiting the qualities for 
quantum computation da and an electrical tunability of 
quasi-flat bands respectively. Recent study on ultracold 
fermionic atoms patterned into finite two dimensional op¬ 
tical lattices of a kagome geometry has added to such 
excitements [l8j |. Advancement in controlled growth of 
artificial, tailor-made lattices having atoms trapped in 
optical lattices with the complications of even a kagome 
structure 0113 has contributed towards increasing in¬ 
terest in these deceptively simple looking systems. 

The non-dispersive character of the energy (E)-wave 
vector (k) curve makes the effective mass of the elec- 
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tron infinite, making the electron having such an energy 
practically immobile in the lattice. The corresponding 
single particle state is sharply localized at a point, or 
in a finite cluster of points in the lattice. Such clusters 
are separated from the neighboring clusters of the lat¬ 
tice by vertices where the amplitude of the wave function 
is zero. The energy eigenvalue corresponding to such a 
construction belongs to a flat band. Fig. [T] displays one 
such case, where the central hexagonal cluster has non¬ 
zero amplitude at the vertices, while the subset of atomic 
sites residing at the vertices of the triangle are zero. This 
figure corresponds to energy E = —2 in a tight binding 
description where the ‘on-site’ potential at every lattice 
point has been set equal to zero, and the nearest neighbor 
hopping integral has been assigned the value unity. Such 
situations can also be realized in the kagome lattice [2l| , 
or other FB systems 221. 
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FIG. 1: (Color online) Schematic view of a part of an el¬ 
ementary star-network illustrating the localized distribution 
of amplitudes within a finite cluster. 


The FB systems have been exemplified so far in terms 
of periodic lattice models. The number of flat bands in 
these systems is always finite. Gomparatively speaking, 
practically nothing is known about the existence of dis¬ 
persionless, flat band modes in hierarchically generated 
artificial lattices, where, because of the inherent scale in- 
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variance, it might be possible to identify an infinity of 
such modes. 

Here we address this issue in terms of a class of Sierpin- 
ski gasket (SPG) lattice^which have been widely stud¬ 
ied over the years [l^- 261 in respect of their electronic 
properties. As one should appreciate, the FB states are 
localized states with finite support in a lattice. Though, 
Anderson localization with weak disorder has been ad¬ 
dressed for tight binding models with flat bands , the 
effect of self similar lattice topology in this context is 
something which we find practically unaddressed. 

This leads to another, and a very important additional 
motivation for the present work, namely, an exact eval¬ 
uation of at least a subsection of localized eigenstates in 
a hierarchical lattice. As such lattices are not periodic, 
one generally has to resort to an exact diagonalization 
of the Hamiltonian to know the energy eigenvalues, and 
because of the fragmented. Cantor set-like energy spec¬ 
trum usually offered by these lattices the eigenvalues 
obtained from one finite generation of the lattice in gen¬ 
eral, slip out of the spectrum on increasing the size of the 
system. A precise knowledge of the energy spectrum in 
the thermodynamic limit thus turns out to be a non triv¬ 
ial problem. Though in recent times the existence of an 
infinity of extended, fully transmitting eigenstates have 
been proven within the framework of real space renor¬ 
malization group (RSRG) schemes [2^, [2^, the exact 
evaluation of localized state-eigenvalues has still lar gely 
remained an unresolved issue, until very recently [30l l3l| 
and, for a special kind of deterministic fractal geometries. 



FIG. 2: (Color online) Portion of an infinite 6 = 3 Sierpinski 
gasket (left) and its renormalized version. The circles col¬ 
ored red and black have coordination numbers four and six 
respectively. The encircled vertices survive after decimation, 
leading to the scaled lattice (right). 


In this communication we focus on the so called b = 3N 
Sierpinski gasket family [s^. We show that, it is possible 
to explicitly ‘construct’ localized eigenstates spanning fi¬ 
nite sized clusters on such lattices. The clusters support¬ 
ing such localized distribution of amplitudes of the wave 
function are separated from neighboring clusters by ver¬ 
tices where the amplitude of the wavefunction is zero at 
the specified energies. As a result, the electron with such 
energies are forced to hop around and remain confined 
in these finite clusters. The minimum size of the clusters 
supporting these localized states increases with the value 


of n. Using an RSRG scheme it is possible to unravel 
an infinity of such localized states and the correspond¬ 
ing spanning clusters for every member of the b = SN 
family. We present explicit results for the 6 = 3 gasket. 
The remaining members can be easily dealt with using 
the same formalism, but with increasingly cumbersome 
algebra. 

In the second phase of our work we examine the effect 
of a perturbation on the FB states in such hierarchical 
networks. The perturbation is chosen in the form of a 
uniform magnetic field piercing each elementary plaque- 
tte in the 6 = 3 gasket. The flux threading each pla- 
quette produces remarkable changes in the energy spec¬ 
trum, generating apparent absolute continuity in the en¬ 
ergy spectrum. Such subbands are populated only by 
extended wavefunctions, as revealed by exhaustive nu¬ 
merical study based on the RSRG scheme. The fractal 
character of the energy spectrum is removed, at least lo¬ 
cally, and the lattice becomes completely transparent to 
an incoming electron when the Fermi level is chosen to 
lie in such continuous subbands. The FB states appar¬ 
ently disappear (at least, they could not be explicitly be 
constructed any more). 

In what follows we describe our findings. In section HIl 
we describe the Hamiltonian and the RSRG scheme. Sec¬ 
tion nil contains a discussion of the construction of the 
localized eigenstates on finite supports. A critical discus¬ 
sion on the flat band status of such states are provided 
in Section |TVl Section IVl highlights the effect of inserting 
a uniform magnetic flux in each elementary plaquette of 
the 6 = 3 gasket, and the two-terminal transport calcula¬ 
tions are discussed in Section IVTl Finally, in section IVTIl 
we draw our conclusions. 


II. THE HAMILTONIAN AND THE RSRG 
SCHEME 


We refer to Fig. [2] where a portion of a 6 = 3 SPG 
network (left) and its renormalized version are shown 
(right). The electron’s hopping is restricted only along 
each edge, and between nearest neighbors only. We dis¬ 
tinguish between two types of sites, viz., one with coordi¬ 
nation number four and colored in red (henceforth called 
C 4 sites) and the other, black circles, having a coordina¬ 
tion number six (Cq sites). The system is described by 
the tight binding Hamiltonian, 


H = ^eic|c* 

(b> 


cjc, -I- h.c. 


( 1 ) 


Ci is the on-site potential at each vertex, and will be as¬ 
signed symbols £4 and £5 for the C 4 and Cq sites respec¬ 
tively. t is the nearest neighbor hopping integral. We 
shall be interested only in the localization induced by 
the topology of the lattice, and hence, in all numerical 
results we shall set £4 = £6 = 0, and t = 1. However, on 
renormalization, the on-site potentials deviate from each 
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other at every stage, and that is why it is necessary to 
distinguish them at the very beginning. 

Using the set of difference equations satisfied by the 
electrons, viz., 

{E - ei)'ipi = ^ tijipj (2) 

3 

we obtain a scaled version of the original lattice, as shown 
in Fig. [5]. An appropriate subset of the and C% sites 
are decimated to yield an identical geometry with the 
renormalized on-site potentials and the hopping integral. 
The ‘surviving’ sites on the renormalized lattice are en¬ 
circled in the left portion of Fig. [21 The RSRG recursion 
relations for the on-site potentials and the hopping inte¬ 
grals are given by, 

^6,n-t-l — “b 3.A Yi 

^4,n-t-l — “b 2.A n 

, 1^2, n.E'^ + 72,+ ^2,n 

tn+1 — -T- 

where, Tn = {ai,nE^ + Pi,nE'^ + ')i,nE -b di,„)/A„, and, 

^l.n ~ 

/3l,n = ~ + 2e4^n + £6,n) 

'Jl.n — ^n(c4,n “b 

C4,n(c4,n “b 2cQ^n')] 

^l.n — ‘2‘^ni‘^^n ~b “b C 4 yjC 6 ,n 

t^(4e4,n + Ce.n)] (4) 

The other symbols used are, 

132,n = tl 

72,n — ^n( 4 ^n ^ 4 ,n ^6,n) 

^ 2 ,n — t^(c4,nC6,n “b ^tn^4,n) 

A„ = [{E - 64,nf - tl][{E - 66 ,,n){E - 64,n “ 2t„) - 6t^] 

( 5 ) 

The ‘stage’ of renormalization is denoted by n > 0, and 
we have chosen to stick to € 4.0 = ee.o = e and to = t at 
the bare scale of length. 

The recursion relations Eq. @ will now be exploited to 
construct and analyze a class of localized and extended 
states in such systems. 

III. CONSTRUCTING THE LOCALIZED 
STATES 

A. The basic scheme 

Let us begin with the 6 = 3 gasket and refer to Fig. (3] 
It is simple to work out that, if we set E = 64 — 2t^ a 
consistent solution of the difference equations (or equiv¬ 
alently, the Schrbdinger equation) is obtained in which 


the amplitude of the wave function ipj is zero at the Cg 
vertices with coordination number six and lying at the 
center of the shaded hexagon, as well as on a subset of 
C 4 vertices. The non-zero amplitudes are designated by 
X (which can be set equal to unity, for example). This is 
illustrated in the left panel of Fig. (3] This latter subset of 
C4 sites sit at the ‘connecting’ points between hexagonal 
plaquettes at different portions of the infinite lattice. By 
making the amplitude vanish at these strategic sites it 
is possible to confine the hopping of the electron around 
the periphery of the shaded hexagons. The eigenstate 
resembles the spirit of a molecular state [ 1 ^ localized 
partially by quantum interference and partially by phys¬ 
ical boundary formed by the sites with zero amplitude. 
The wave function corresponding to this energy eigen¬ 
value does not have any overlap with the wave function 
described on the next neighboring hexagons. The local¬ 
ization area is the shaded hexagonal plaquette confining 
the electron’s motion for E = 64 — 2t along the ‘periph¬ 
ery’, and is distributed throughout the infinite geometry. 

B. Exploiting the self similarity 

Once such a construction of localized state-amplitudes 
is made possible in the bare length scale, the self sim¬ 
ilar structure of the SPG lattices can be exploited to 
construct wave functions, and to ‘work out’ the corre¬ 
sponding areas of localization at various scales of length. 
The eigenvalues corresponding to such cases are easily 
extracted from the solutions of the equation 

E = e4,n - 2tn (6) 

at any desired n-th stage of RSRG. For every real root 
of this equation, a construction similar to that in the 
original lattice can be made in respect of the values of the 
amplitudes of the wave function at different C4 and Cg 
vertices on the scaled version of the lattice. One then has 
to work in the reverse direction to lay out the amplitude- 
distribution on the original lattice at the bare scale of 
length by simultaneously solving the difference equations 
Eq. ([2]) for the respective energy eigenvalue. The scheme 
is illustrated in Fig. EKb). With 64 = cq = 0 and t = 1, 
we get E = — 2 at the original length scale, while the 
same parameters yield localized states at E = —2 and 
E = 1 ± \p2 when we set n = 1, that is, when we solve 
the equation E = £4 4 — 2ti. Let us now explain the 
construction of localized states for the new roots, viz., 
A = I ± V 2 . 

The hexagon in Fig. EKb) whose sides are highlighted 
by thick green line is the ‘smallest’ hexagon on a one step 
renormalized gasket. We choose the distribution of am¬ 
plitudes at E = £ 4.1 — 2 ti around this hexagon exactly in 
the same way as we did for A = —2 in the bare length 
scale (Fig. jSKa)). Such a construction, as we now appre¬ 
ciate, can be made consistently throughout the infinite 
system. The next job is to work out, in the reverse direc¬ 
tion, the amplitudes of the wave function at all the ver- 
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(a) 



FIG. 3: (Color online) (a) Distribution of a fiat-band state amplitudes on a 6 = 3 SPG network at its bare scale of length, 
and (b) the distribution on a one step renormalized gasket. The localized state eigenvalues are extracted from the equations 
E = 64,0 — 2to and i? = 64,1 — 2t\ respectively. 


tices of the original gasket with the prefixed amplitudes 
(with values ±x at the edges of the green highlighted 
hexagon and zero at its center) for energy eigenvalues 
E = l±V2. 

The amplitudes are denoted by the letters a-f in 
Fig. [31(b). We have shown explicitly the smaller hexagons 
(blue highlighted) that are embedded in the larger one 
(green highlighted). It is possible to obtain a consis¬ 
tent set of values for the amplitudes at all vertices of the 
smaller (blue) hexagons. For example, with x = 1, set 
arbitrarily, the distribution of amplitudes for E = 1 + \/2 
(one solution of if = £ 4,1 — 2 <i with e = 0 and t = 1) is 
given by, a = —/ = (v^ — l)/2, b = —e = l/-\/2 and 
c = —d = 1/2. The amplitude at the center of a small 
hexagon is zero for this particular distribution. An iden¬ 
tical pattern of distribution, but with a different set of 
values for a, 6, ... /, can be easily obtained if one sets 
A = 1 - ^2. 

Needless to say, the scheme can be extrapolated to 
unravel other localized states by solving the equation 
E = £ 4,71 — 2tn at any desired level. The roots increase 
in number, and the set of solutions obtained at the n-th 
stage of renormalization contains all the roots obtained 
from the (n — l)-th stage. In every case one has to iden¬ 
tify the smallest hexagon in that scale of length, set the 
amplitude of the wave function equal to zero at its center. 
The sites decorating the periphery of the basic hexagon 
then should be assigned amplitudes equal to Ex as we 
did in Fig. |3 Such a basic hexagon at the renormalized 
scale must be separated from its neighboring hexagons, a 
condition that is easily satisfied by setting the amplitude 
at the ‘connecting’ sites equal to zero. The final task is 
to map the distribution back onto the original lattice and 
to solve the set of equations Eq. ([2]) to work out the am¬ 
plitudes at the inner sites of the lattice at its bare scale 


of length. 


C. The ‘staggered’ character of localization 

An interesting qualitative distinction between the en¬ 
ergy eigenvalues corresponding to these localized states 
can be unravelled if we look at the evolution of the hop¬ 
ping integral under RSRG steps. For any energy eigen¬ 
value obtained from the n-th stage, and by solving the 
equation E = £ 4 ,„ — 2t„, the hopping integral remains 
finite upto n-th RSRG steps. This means that one has 
non-zero overlap of the wave function between nearest 
neighboring sites at that length scale. When mapped 
back onto the original lattice it amounts to a non-zero 
overlap of the wave function between sites much beyond 
the nearest neighbors. However, the hopping integral 
starts decaying after the n-th RSRG step implying that 
the wave function looks ‘extended’ only locally, and in re¬ 
ality such wave functions are localized on a finite cluster 
of the infinite lattice. The overlap of the envelopes of the 
wave function for a localized state-energy extracted at 
the n-th stage weakens beyond a certain scale of length. 
Such finite sized clusters span larger and larger areas in 
real space as one extracts the energy eigenvalue by solv¬ 
ing the above equation at deeper and deeper scales of 
length.Thus simply by ‘choosing’ the value of stage of 
renormalization ’n’ one can delay, in space, the localiz¬ 
ing character of the wave function. In a recent study [s^ 
we encountered such staggered localization in a special 
class of fractal lattices. 

Before we end this section, it is pertinent to point out 
that the size of the spanning clusters even at the bare 
length scale depends on the value of n in the class of 
b = 3N SPG fractals. Fig. H] depicts three other cases. 


























(a) 


(b) 


(c) 


FIG. 4: (Color online) Distribntion of amplitudes for a non-dispersive eigenstate at E = €4 — 2t on b = 6 , b — 9 and b = 12 
Sierpinski gaskets. With increasing scale factor, the amplitudes assume non-trivial distribution over larger and larger clusters 
which ultimately get decoupled from the rest of the lattice, courtesy, the zero amplitudes at the connecting vertices. 


viz., b = 6 , b = 9 and b = 12. We stick to the same value 
of energy, viz., E = 64 — 2t. The localization areas even 
at the bare length scale are already larger than that in 
the b = 3 case. In each case, the electron’s hopping is 
confined around the peripheries of clusters of finite size, 
shown by the shaded areas. Such finite sized clusters 
grow bigger and bigger as one goes from b = 6 to b = 12, 
and beyond. The clusters, however big, are ultimately 
separated from each other by a subset of C 4 sites with 
ipj = 0 , and are distributed throughout the infinite ge¬ 
ometry. 


IV. THE FLAT BAND ANALOGUES 

A pertinent question is, whether these infinite number 
of localized eigenstates are non-dispersive in character, 
in the spirit of the fiat band states observed in a class of 
periodic lattices as mentioned in the introduction. The 
absence of translational invariance in the SPG class of 
lattices prevents a derivation of the usual energy-wave 
vector dispersion relation. At the same time, the fact 
that the localized states discussed so far have been ‘con¬ 
structed’ so as to display a kind of self localization^ where 
the amplitudes are distributed over a finite cluster of lat¬ 
tice points, effectively ‘decoupled’ from the rest of the 
system. This inspires us to undertake a couple of indi¬ 
rect, but effective ways of discerning the actual character 
of these states. 


A. Dispersion relation for a periodic array of 
fractal approximants 

First, we construct a periodic array of 6 = 3 SPG net¬ 
works. Each ‘unit cell’ consists of a finite SPG at its n-th 
generation. One such array comprising the first genera¬ 
tion 6 = 3 SPG network is shown in Fig.[S]for illustration. 


As is evident, such a construction requires the presence 
of a top vertex with coordination number two (colored 
cyan in Fig. 0. It is important to appreciate that for 
such finite systems the distribution of amplitudes does 
not follow the method prescribed for Fig. 0 which was 
typical of an infinite system where the ‘end points’ were 
not ‘visible’. 

As we gradually increase the value of n, the unit cell 
itself turns out to be a reasonably good approximant of 
the infinite system, and becomes the true representative 
in the limit n —^ 00 . For a periodic array of n-th gener¬ 
ation SPG networks one can renormalize each ‘unit cell’ 
to n-times to convert the array into an array of trian¬ 
gles with the potentials at the vertices equal to C 2 ,n and 
C 4 ^n respectively. The nearest neighbor hopping integral 
along the edges of such effective triangles are tn- This ar¬ 
ray of triangles is then converted into an array of single 
‘effective’ atomic sites (encircled by dotted lines in Fig. [5]) 
with energy dependent potential, and inter-connected by 
an energy dependent hopping integral. The effective, en¬ 
ergy dependent on-site potential and the nearest neigh¬ 
bor hopping integral in the linear chain of renormalized 
atoms are given by. 


2 tl 


e = e4,i 


t = 


E — £2,; 


E — £2 . 


(7) 


In this case, we need an additional renormalization of the 
top vertex with coordination number two, viz., 

e2,n+l = ^2,n + iFn 

The definition of Fn is already given in earlier section 
(Eq.(|3]) and Eq.®). The E-k dispersion relation for such 
a periodic chain is then easily obtained through the equa¬ 
tion, 

E = e + 2icos ka 


( 8 ) 
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FIG. 5: (Color online) A periodic array of finite generation b = 3 SPG networks. With increasing generation each unit cell will 
be representative of an infinite SPG network. 


In Fig. [ 6 ] we depict the dispersion relation when the 



k k 




FIG. 6: (Color online) The dispersion relation for the array 
depicted in Fig. [5] The flat, dispersionless modes are ap¬ 
parent. The flat band at E — —2, the root of the equation 
E = 64 — 21 with £4 = 66 = 0 and t — 1 has already been shown 
to exist, by explicit construction, for an infinite system. 

unit cell consists of (a) the first generation, (b) the second 
generation, (c) the third generation and (d) the fourth 
generation 6 = 3 SPG networks. We have set €4 = €q = 0 
and t = 1. It is clearly seen form Fig. EKc) and (d) 
that flat, non-dispersive bands appear at E = —2 and at 
E = 1± 1/2 - the three roots obtained from the equation 
E = — 2tn at the bare length scale and from the one 

step renormalized lattice, as discussed in the last section. 
Even with a modest unit cell structure for the triangular 
array, the self (staggered) localized states corresponding 
to the truly infinite system are showing up. Thus, we 
feel tempted to identify these self-localized states as flat 
band analogues of similar states in quasi one dimensional 
systems [lo|-[i3]- However, we are now in a position to 
appreciate that such states will be infinite in number as 
the size of the system approaches the thermodynamic 
limit. 


B. A study of the density of states 

As a second check, we evaluate the local density of 
states (LDOS) at the C 4 and Cq sites. The results are 


shown in Fig. [TKa) and (b) respectively. This is done 
following the standard procedure, and using the RSRG 
recursion relations Eq. ©, which finally yield the lo¬ 
cal Green’s function (LGF) at the C 4 and Cq sites as, 
~ £ 4 )”^ ~ respectively. 

The symbol e* indicates the fixed point value of the re¬ 
spective parameter. Finally, the LDOS at the desired site 
is obtained as, pj = (—I/tt) lim,)_,.o Im GQQ{E + ir]), with 
j = 4 or 6 depending on the coordination number. 

The LDOS (in the absence of an external magnetic 
field) exhibits a highly fragmented distribution as seen 
in Fig.[7Ka) and (b). This is expected. The fine structure 
is related to the distribution of spectral weight [sj, and 
is controlled by the lattice topology. One should observe 
that the j-th site {C 4 or Cq) having = 0 can not con¬ 
tribute to the LDOS at that particular site. Nevertheless, 
non-trivial contribution comes from the neighboring sites 
yielding a finite value of the density of states at a C 4 site. 

It is important to appreciate the role of the numeri¬ 
cal value of the imaginary part rj added to the energy. 
77 determines the width of the LDOS around an energy. 
We have carefully examined the growth of the LDOS for 
E = —2, 1 ± \f2 and at several other localized flat band 
states as obtained by explicit construction at bare length 
scale and from a one-step renormalized lattice. The result 
is tabulated in Table U for the above three energy eigen¬ 
values. It is explicitly seen that, for each of these energy 
eigenvalues, with a sequential decrease in the value of 77 , 
the LDOS at the C 4 site, viz., p 4 increases exponentially, 
indicating that in the limit 77 —> 0, the LDOS diverges 
at these special energy values. As we know, the average 
density of states (ADOS) is related to the group velocity 
Vg of the electron by the relation 

Pav oc j v~^df‘k (9) 

where, Vg oc dE/dk. Therefore, a diverging LDOS, lead¬ 
ing to a diverging ADOS, is indicative of the zero (or ex¬ 
tremely low) group velocity and hence, a non-dispersive 
character of the state. The energy eigenvalues, obtained 
by construction, and corresponding to the self-localized 
eigenfunctions in the present example can thus be con¬ 
sidered as the so called flat band states. We can thus 
unravel infinite number of such non-dispersive modes in 
an infinite b = 3N SPG network. 
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E 

V 

P4(E) 

pa{E) 


10 “^ 

3.28 

0.03 


10 “® 

15.66 

0.07 



74.38 

1.5 X 10 "® 

-2 

10 "® 

355 

3.2 X 10 "'‘ 


10 “® 

1703.2 

6.8 X 10 "® 


10 -^ 

8701 

1.34 X 10 “® 


10 “® 

40419 

2 X 10 "® 


10 “^ 

1.68 

1.06 


10 “® 

5.2 

0.5 


lO"'^ 

22.5 

0.11 

l-y/2 

10 “® 

106.4 

2 X 10 "^ 


10 "® 

506.24 

5.4 X 10 "® 


10 -" 

2402.8 

1.13 X 10 “® 


10 "® 

11454 

2.4 X 10 "^ 


10 "^ 

1 

7.5 X 10 "® 


10 “® 

4.41 

1.75 X 10 "® 


10 "“ 

20.91 

3.8 X 10 "® 

1 + ^/2 

10 “® 

100 

8 X 10 "'‘ 


10 "® 

473 

1.7 X 10 "'‘ 


10 -" 

2249 

3.6 X 10 "® 


10 “® 

10747 

7.4 X 10 "® 


TABLE I: The values of the local densities of states Pa[E) 
and pe,{E) at special energy eigenvalues E = —2, and 1 ± y/2 
as obtained from the original lattice and its one-step renor¬ 
malized version corresponding to localized flat-band states. 
We have set the on-site potential e = 0 at each site, and the 
energy is measured in unit of the hopping integral t. 


V. EFFECT OF A UNIFORM MAGNETIC 
FIELD 


As is well known, the FBS are unstable against pertur¬ 
bation. We examine the general character of the energy 
spectrum of a 6 = 3 SPG network when each elemen¬ 
tary plaquette is threaded by a uniform magnetic flux $. 
The presence of a magnetic field breaks the time rever¬ 
sal symmetry of the electron-hopping along a bond, and 
this is taken care of by introducing a Peierls’ phase in 
the hopping integral, viz., t —t exp (*27r$/3<i>o), where 
$0 = hcje is the fundamental flux quantum. 

The RSRG method is used again to evaluate the 
LDOS, which is shown in Fig. [TJ^c) and (d) for the 
and Ce sites respectively for $ = $o/4. We draw the at¬ 
tention of the reader to the absolutely continuous shaded 
portions in the LDOS. This is a generic character for all 
values of the piercing flux. We have carefully scanned 
the continuous subbands in arbitrarily small energy in¬ 
tervals, and observed the flow of the hopping integral 
under the RSRG iterations. Extensive numerical search 
reveals that for any energy belonging to the continua, 
the hopping integral keeps on oscillating for an indefi¬ 
nite number of iteration loops. This implies that, for 
any energy forming the continua there is non-zero over¬ 


lap between the envelopes of the wave function at various 
scales of length, and that, the subbands are densely pop¬ 
ulated by extended wave functions. In addition, the en¬ 
ergy E = 0 lying at the center of the spectrum exhibits 
a one cycle fixed point in the parameter space, that is, 
(^4,n j ^ 6,717 \ tn \) — (^4,n-t-l 7 ^ 6 ,n-t-l 7 |^n-t-l | ) 7 7Z ^ 0. 
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FIG. 7: (Color online) Local density of states at the (a) G 4 
and (b) Cq sites of an infinite b = 3 SPG network in the 
absence of a magnetic field, (c) Local density of states at a Ca 
site when each elementary plaquette is pierced by a magnetic 
flux $ = 4>o/4, and (d) the same at a Ce site with ^ — ^0 /A. 
We have set £4 = eg = 0 and t = 1. The imaginary part added 
to the energy is rj — 10 ~^. 


The existence of the subbands of extended states does 
not rule out the other localized states, though the non- 
dispersive, flat bands apparently do not exist (or, at least, 
become undetectable). As, for an electron circulating 
around a closed loop, the mag netic flux plays a role equiv¬ 
alent to the wave vector 3^, we illustrate the E-^ dia¬ 
gram (Fig. [8]) as an analogue of the dispersion relation 
in Fig. [51 The energy eigenvalues of the fractal network 
are seen to form mini bands as a function of the flux <I> 
with period $ 0 . There are multiple, inter-twined band 
crossings, and quite a dense distribution of eigenvalues, 
forming almost continuous quasi-flat E-^ ‘bands’. The 
pattern is flux periodic with a period equal to one flux 
quantum, and close observation of the E-^ ‘band’ sub¬ 
sections reveals formation of interesting variants of the 
Hofstadter butterflies [s^. The spectral landscape is a 
quantum fractal, and encoding the gaps with appropriate 
topological quantum numbers remains an open problem 
for such deterministic fractals. 


VI. TWO-TERMINAL TRANSMISSION 
CHARACTERISTICS 

For the sake of completeness of the above discussion, 
we have computed the two-terminal transmission charac¬ 
teristics of a finite system. The procedure is standard and 
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FIG. 8 : (Color online) Distribution of energy eigenvalues 
against the magnetic flux $ (in unit of <I>o = hc/e) in an 
infinite 6 = 3 SPG network. 


is often used to evaluate the transmission coefficient of 
such hierarchical fractal structures . The basic idea is 
to clamp the system in between two semi-infinite periodic 
leads, the so called ‘source’ and the ‘drain’. The on-site 
potential at the atomic sites in the leads is Ul ■ The hop¬ 
ping integral is 7 l in between the nearest neighboring 
atomic sites in the leads. The finite sized network sand¬ 
wiched in between the two ordered leads is then succes¬ 
sively renormalized to reduce it to an effective diatomic 
molecule (dimer) [s^ . The transmission coefficient of the 
lead-network-lead system then is given by a well-known 
formula [s^, 


with, 


4 sin^ ka 
^ \A\^ + \B\^ 

A = [{Mi 2 — M 21 ) + {Mil ~ M 22 ) cos ka] 
and B = [{Mu + M 22 ) sin ka] 


( 10 ) 


where, refer to the dimer-matrix elements, written 
appropriately in terms of the on-site potentials of the 
final renormalized left (L) and right (R) atoms at the ex¬ 
tremities of the finite SPG network and the renormalized 
hopping between them, cos ka = {E — Ul )and a 
is the lattice constant in the leads which is set equal to 
unity throughout the calculation. 

It is apparent from Fig. |9]that, the low transmitting 
character of the gasket in the absence of an external mag¬ 
netic field (Fig. UKa)) is enriched by high transmitting 
windows in the presence of a magnetic flux (Fig. EKb)). 
This corroborates the argument about the existence of 
absolutely continuous subbands populated by extended 
eigenfunctions. 

The transmission coefficient, as expected, is flux peri¬ 
odic, with period equal to one flux quantum. The fine 
structures of the Aharonov-Bohm oscillations in trans¬ 
port are of course, sensitive to the chosen energy eigen- 



FIG. 9: (Color online) Transmission coefficient of a 4-th gen¬ 
eration 6 = 3 SPG network with (a) $ = 0, and (b) <I> = "Fo/d. 
We have chosen 64 = ee = 0, t = 1, Ul = 0 and Tl =?>■ 


value. However, to save space, we refrain from showing 
these figures. 


VII. CONCLUSION 


In conclusion, we have unravelled an infinity of flat 
band-like, dispersionless eigenstates in the b = 3N class 
of Sierpinski gasket fractals. The states are localized in 
clusters of finite size, which are effectively separated from 
similar clusters in the lattice by atomic sites where the 
amplitude of the wave function is zero. The topology 
of the lattice confines the dynamics of the electron in 
such local clusters for these special sets of energy eigen¬ 
values. The states have been obtained by explicit con¬ 
struction and the inherent self similarity of the lattice is 
exploited through a real space rescaling of the system. 
The dispersionless character of the flat band states are 
conhrmed by obtaining the E-k dispersion curves of a pe¬ 
riodic array of the approximants of the true gasket, and 
also from the diverging character of the local density of 
states at those special energy eigenvalues. The two termi¬ 
nal transport and the magnetic field induced absolutely 
continuous subbands have been discussed in details. The 
intricate mixture of complexity, frustration and a deter¬ 
ministic geometry may throw up new challenges in the 
quantum Hall physics in such fractal geometries. 
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